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The production of photons in heavy ion collisions is dicussed. 
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It has long been thought that the production of photons could be used 
to detect the formation of the quark-gluon plasma (QGP) in ultrarelativistic 
heavy ion collisions. In fact, it is believed by some that the present WA98 
data lH] on inclusive photon production are evidence for QGP formation. At 
CERN, the available center of mass energy of the nucleon-nucleon system 
is y/s ~ 17 GeV. Since the rate of photon emission from the QGP phase is 
expected to increase with energy, relatively to other production mechanisms, 
it should become more and more relevant as one goes to RHIC (-^/s ~ .2 TeV) 
and LHC ~ 5.5 TeV). 

The production of photons in heavy ion collisions is rather complex and 
one roughly distinguishes four mechanisms: 

— 1) the photon is produced in the hard interaction of two partons 
in the incoming nuclei similarly to the well-known Chromodynamics (QCD) 
processes (QCD Compton, annihilation, bremsstrahlung) in nucleon-nucleon 
collisions. The rate is calculable in perturbative QCD and falls off at large 
transverse momentum, Pj,, as a power law; 

— 2) in the collision of two nuclei the density of secondary hadrons is so 
high that quarks and gluons become unconfined and a bubble of hot quark- 
gluon plasma is formed: at LHC the temperature is expected to be of the 
order of 1 GeV. Photons are emitted in the collisions of quarks and gluons 
in the plasma with an energy spectrum which is exponentially damped but 
which should extend up to several GeV; 

— 3) the QGP bubble expands and cools until a temperature of 150 to 
200 MeV is reached and a hadronic phase appears. As they collide the hot 
hadronic resonances (vr'^, p, uj) emit photons until the freeze-out temper- 
ature is reached. The typical energy of such photons ranges from several 
hundred MeV to several GeV; 

— 4) photons are decay products of resonances (vr'', r], etc.) emerging at 
the end of the thermal evolution (in which case their energy is of the order 
of a few 100 MeV). vr'^ and r] resonances can also be produced at large Pj, 
in hard parton scattering at the beginning of the collision (their energy is 
of the order of a several GeV). These photons together with those in 1) are 
a background to thermal photons produced in 2)-3). 

A quantitative application of the above to the WA98 results leads to the 
conclusion that it is ... impossible to conclude whether a QGP is formed! 
Indeed the uncertainties on each mechanisms are such that the data can be 
accomodated by production of type 1) only (the "German" school 0) or a 
mixture of type 1) and type 2-3) (the "Indian" school and the "Finnish" 
school IQ]). The problem is that neither data nor theory tell us how to 
reliably normalise, in the Pj, range of a few GeV, the rates from the various 
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mechanisms. In the fohowing we consider only inclusive photons of energy in 
the GeV range where thermal production (the "signal") competes with the 
other mechanisms and we do not deal with photon-hadron or photon-photon 
correlations which can also be used to characterise the formation of QGP. 
We concentrate on the discussion of uncertainties in the partonic production 
on the one end, and on the rate of production in the quark-gluon plasma on 
the other hand. A review on uncertainties associated to the hydrodynamic 
evolution of the hot matter can be found in recent works, in particular [||, ^. 



1 Partonic production of photons 

Since this mechanism is the same as that in proton-proton collisions, a lot 
of what is known for this process should be relevant for heavy ions. The 
production cross section of direct photon in perturbative QCD is well known 
and has the usual factorisable form 0: 



dp^dy 



/dz 
dxadxb^F^^^ixa, M)F^^^ (xfe, M)D^/k{z, M^) 



da 



dpk^dyk 

where the functions F are the parton densities and a^^ ^ is the hard cross 
section between partons i and j, in the nuclei A and B, respectively to pro- 
duce parton k. The function D^i^ is the fragmentation function of parton 
k into a photon. It reduces to a (5(1 — z) function in case the photon is pro- 
duced directly {k = 7) in the hard collision: in this case one refers to direct 
production of the photon, otherwise one refers to bremsstrahlung produc- 
tion. The calculations have been carried out up to next-to-leading order in 
QCD and they suffer from the usual rather large ambiguities in the choice of 
the renormalisation scale /i, the factorisation scale M and the fragmentation 
scale Mp Resumming "threshold corrections" which is appropriate for 
very large p^, reduces the dependence on scales Q. Extensive phenomeno- 
logical studies have been carried out for proton/(anti-)proton scattering for 
^/s from 20 GeV to 1.8 TeV and the situation is very confused. Concerning 
data on inclusive photon production it is impossible to accomodate all data 
within with the same set of parameters: the E706 data (on Beryllium at 31.6 
and 38.8 GeV) are at least a factor 2 to 3 above the other data which extend 
from 20 to 63 GeV, and they do not have the same Pj. dependence §. There 

^We do not discuss here the impact parameter dependence of the coUision. 



2 



is therefore little chance that present experiments can be used to normalise 
predictions for future heavy ion collisions. 

Several further problems are specific to heavy ion collisions. They fall in 
two categories. 



1.1 Further theoretical uncertainties 

The kinematical domain of interest in heavy ion collisions is not exactly the 
same as the domain where theory has been confronted to experiments. For 
WA98 data, both the energy and the relevant Pj, range are rather low: is 
perturbation theory reliable there? The present data in the same energy 
range do not help either (E629 and NA3 Q) as they have some relative 
normalisation problems. For RHIC and LHC the interesting values of p^, to 
distinguish mechanism 1) from 2-3) above, are a few GeV, i.e. very low p^,. 
This introduces a ratio of large scales (p^ / s) which will become very small 
and the standard NLO calculations are probably not reliable then. This is 
the domain where "recoil" resummation is important but further studies are 
needed since the present results are rather dependent on non-perturbative 



parameters |10|. 



In the small p^, region, at high energies, the production of photons by 
bremsstrahlung becomes overwhelming (for LHC, around 70% at ~ 5 
GeV/c decreasing to 40% at p^ ~ 50 GeV/c). However, the fragmentation 



functions of a quark or a gluon into a photon |11] are not reliably known 
in the relevant kinematical domain. Preliminary studies show a variation 
of a factor 4 in the predicted cross section (at small p^) depending on the 
two fragmentation sets used, both of which beeing equally good fits to lower 
energy data [p^]. 



1.2 Nuclear effects 

Three effects can be considered: shadowing, "primordial A:^" and parton 
energy loss. 

If shadowing is not a problem for fixed target experiments, such as WA98, 
because they probe rather large x values of the partons in the nuclei, it 
is however predicted that it will reduce, at RHIC and LHC, the rate of 
production of photons below = 50 GeV since the useful x range (namely 
X < .1) is well within the domain where shadowing occurs: at LHC the 
suppression can be of the order of 30% in the small p^, region |13|. Gluon 



initiated processes play an important role, however not much is known about 
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gluon shadowing and to make predictions one has to rely on models [^] 
which introduce further uncertainties in the predictions. 

Nucleons in nuclei are expected to have some transverse momentum with 
respect to the beam axis. Assuming that this transverse momentum is com- 
municated to the partons which scatter in type 1) processes, one has to 
modify eq. (1) to take it into account. The effect is incorporated in leading 
logarithm calculations in a simple way. In the usual implementation p5|] , 
the main effect is to modify the partonic cross section a which behaves as 
in eq. (1) to — fc^' — A;*"') where the fc^"' are the intrinsic trans- 

verse momenta of the incoming partons. Such a modification considerably 
increases the size of a and consequently the photon production rate spe- 
cially if is comparable to fe^"'' . In fact, a cut-off is introduced to prevent 
the vanishing of the denominator. To a large extent, the primordial kj, effect 
is controled by the choice of this arbitrary cut-off. Clearly, efforts should be 
made to improve the method and provide a basis for reliable extrapolations. 

As mentioned above, at RHIC and LHC, bremsstrahlung production of 
photons dominates the rate at low Pj, . But the final state parton looses en- 
ergy jl^ in the hot medium "before" emitting the photon, thereby affecting 
the energy spectrum of the photon and reducing its rate. In a recent model 
calculation it was found that the reduction could be, very large, as much as 
70% for LHC at = 3 GeV/c, becoming less drastic as Pj, increases |13]. 
Further studies are necessary to include the latest theoretical results and 
reduce the model dependence. 



2 Photon production in the quark-gluon plasma 

We assume the plasma in thermal equilibrium (temperature T) with van- 
ishing chemical potential. In principle, we could calculate the rate of pro- 
duction in the plasma using an expression similar to eq. (1), substituting to 
the nucleus structure fonctions the Fermi-Dirac distribution for the quark 
and the Bose-Einstein distribution for the gluon, not forgetting the appro- 
priate suppression or enhancement factors for partons in the final state. The 
matrix elements for the various processes (QCD Compton, annihilation and 
bremsstrahlung production) remain the same as before. A crucial difference 
with eq. (1) is the phase space: there, the beam defines a longitudinal direc- 
tion and requiring the photon to be produced at large transverse momentum 
guarantees that the internal lines in the hard matrix elements are finite (the 
denominators in the matrix elements never vanish) and the production rate 
remains finite. In the plasma, because quark and gluon distributions are 
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isotropic, there is no longitudinal direction with respect to which a mini- 
mum transverse momentum for the photon can be imposed, and a naive use 
of eq. (1) leads to a divergent result. The solution to this difficulty comes 
from the fact that the partons in the plasma are not free but are in inter- 
action with other quarks and gluons in the medium: their properties are 
affected over long distances or equivalently when they have soft momenta. 
The theoretical framework to incorporate these effects is the Hard Thermal 
Loop (HTL) effective theory 0. 

Following the HTL approach one distinguishes two scales: the "hard" 
scale, typically of order T or larger (the energy of quarks and gluons in the 
plasma) and the "soft" scale of order gT where g, the strong coupling, is 
assumed to be small. Collective effects in the plasma modify the physics 
at scale gT i.e. over long distances of 0{l/gT). These effects lead to a 
modification of the propagators and vertices of the theory and one is led 
to introduce effective (re-summed) propagators and vertices. This is easily 
illustrated with the example of the fermion propagator, S{P), which in the 
"bare" theory is simply 1/p (we neglect spin complications and make only a 
dimensional analysis). The thermal contribution to the one loop correction 
S(p) is found to be T,{p) ~ g'^T'^/p which is of the same order as the inverse 
propagator when p is of order gT. The re-summed propagator *S{P) = 
l/{p — S(p)) is then deeply modified for momenta of 0{gT). Likewise, the 
gluon propagator and vertices are modified by hard thermal loops when 
the external momenta are soft. A practical consequence of these thermal 
effects is to provide physical infrared cut-offs (effective masses) to quarks 
and gluons which regularise the divergences mentioned above in the naive 
approach. It also leads to new processes which increase the rate of photon 
production One can construct an effective Lagrangian [19| in terms of 
effective propagators and vertices and calculate observables in perturbation 
theory. 

The rate of production, per unit time and volume, of a real photon of 
momentum (i?, p) is calculated via the formula: 

^ = -pi)3«.(^)i-nV(£.rt. (2) 

where ^^{E,p) is the retarded photon polarisation tensor. The pre-factor 
ng{E) provides the expected exponential damping exp{—E/T) when E ^ 
T. A loop expansion of H^ is constructed with effective propagators and 
vertices. 

At one-loop, the photon production rate occurs via the usual QCD 
Compton and annihilation processes the quarks having now an effective 
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Figure 1: Physical processes occuring at two-loops. I: bremsstrahlung from 
an antiquark. II: qq annihilation with scattering. Ill: bremsstrahlung from 
a quark. 



mass. The result is finite and can be written as |20|: 

mi 



Im {E,p)^ e'g'T' {ln{^) + c) (3) 



,2 _ , ,,2rr2 



where m^ ~ g T is related to the thermal mass of the quark which acts as 
the "soft" cut-off . This leads to a "large" logarithmic term of type ln(l/g() 
dominating over a "constant term" C. The two- loop diagrams contribute 
the processes shown in the figure. The result for hard photons is 



luilf{E,p) r^e^g^T"^ (4) 

brcms 



lmU''{E,p) ^ g^TE (5) 

annil 

These two-loop contributions have the same order as the one-loop ones be- 



cause of strong "collinear" singularities |21] which manisfest themselves as 
a factor of the form T^/m^ ~ l/ff^: this factor compensates the extra g"^ 
factor in the numerator. Another interesting result of the calculation is the 
importance of process II in the figure: it describes qq annihilation with scat- 
tering (a "3 to 2" process crossed from bremsstrahlung) which grows with 
the energy of the photon and dominates over the other contributions when 
E /T ^ 1 . Phenomenological applications of these results have been carried 
out and the two-loop processes have been included in hydrodynamic evolu- 
tion codes to predict the rate of real photon production. It is found that 
the two-loop processes (especially the annihilation with scattering) lead to 
a large increase of the rate ^, ^, ^ . 

Since the one-loop and two-loop results are of the same order it is rea- 
sonable to worry about the convergence of the perturbative expansion in 
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the effective tfieory. Tfie enliancement mecfianism operative at two-loop 
could also be at work at the multi-loop level especially in ladder diagrams. 
It was indeed found that the higher order ladder diagrams contribute to 



the same order as the two- loop ones [22|. Another effect which can mod- 



ify the collinear enhancement mechanism is related to the fermion damping 
rate: it can be incorporated by including a width on the fermion propagator 
which takes into account the interaction length of the fermion in the plasma. 
This is technically a higher order effect but it nevertheless modifies the rate 



of production at the leading order |22]. It is interpreted as the Landau 



Pomeranchuck-Migdal effect [ p^]. Both these higher order effects have been 
combined in a consistent way ||2J] and it is found a reduction of the 2-loop 
photon production rate. From the phenomenological point of view the re- 
duction is less that 30% in the Pj, range of interest. All these calculations 
are done for a plasma in chemical equilibrium which may not be a realistic 



approximation: the effect of fugacities has been considered in |25, 26 1: non 
linear effects occur since decreasing the fugacities will reduce the rate of 
photon production but, because the effective fermion mass is also reduced, 



the collinear enhancement becomes stronger [26| 



In conclusion, the qualitative picture of the evolution of heavy ion col- 
lisions seems rather well understood but a lot of progress has to be made 
for the model to become really predictive. Only a few of the open problems 
have been discussed here and there are many more such as the problems of 
thermalisation and chemical equilibration which relate to the initial state 
conditions before the hydrodynamic evolution takes place. 
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